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ABSTRACT 

We describe and test a method to quantitatively classify clusters of galaxies 
according to their projected morphologies. This method will be subsequently used to 
place constraints on cosmological parameters (D and the power spectrum of primordial 
fluctuations on scales at or slightly smaller than that of clusters) and to test theories 
of cluster formation. We specifically address structure that is easily discernible in 
projection and dynamically important to the cluster. The method is derived from the 
two-dimensional multipole expansion of the projected gravitational potential and yields 
dimensionless power ratios as morphological statistics. If the projected mass profile is 
used to characterize the cluster morphology, the power ratios are directly related to the 
cluster potential. However, since detailed mass profiles currently exist for only a few 
clusters, we use the X-ray-emitting gas as an alternative tracer of cluster morphology. 

In this case, the relation of the power ratios to the potential is qualitatively preserved. 

We demonstrate the feasibility of the method by analyzing simulated observations 
of simple models of X-ray clusters using the instrument parameters of the ROSAT 
PSPC. For illustrative purposes, we apply the method to ROSAT PSPC images of 
A85, A5I4, AI750, and A2029. These clusters, which differ substantially in their X-ray 
morphologies, are easily distinguished by their respective power ratios. We discuss the 
suitability of this method to address the connection between cluster morphology and 
cosmology and to assess whether an individual cluster is sufficiently relaxed for analysis 
of its intrinsic shape using hydrostatic methods. Approximately 50 X-ray observations 
of Abell clusters with the PSPC will be amenable to morphological analysis using the 
method of this paper. 
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1. Introduction 


Clusters of galaxies, the largest bound objects in the universe, are especially useful laboratories 
for probing the underlying cosmology (for reviews see Oegerle, Fitchett, & Danly 1990; Durret 
et al. 1995). The dynamical state of clusters in particular has received much recent attention 
because of its implications for estimations of the detailed mass distributions (e.g., Fitchett 1990) 
and intrinsic shapes (Buote & Tsai 1995a) of clusters, and it may provide a powerful constraint 
on the cosmological density parameter (Richstone, Loeb, & Turner 1992, hereafter RLT; Evrard 
et al. 1993). These studies emphasize that the dynamical state of a cluster is qualitatively related 
to the degree of substructure present. Hence there is a need for a scheme that quantifies the 
morphology of clusters in relation to how much “dynamically relevant” substructure they possess. 

What constitutes “dynamically relevant” substructure is at present poorly defined. For 
cosmological purposes RLT make the provisional suggestion that the relevant substructure must 
have a density contrast compared to the mean cluster density in the range 2-10, contain at least 
20% of the cluster mass, and lie within a projected radius of 1.5 Mpc. Denser structures would 
also be important for studies of the underlying mass distribution and intrinsic shape. 

Present techniques to measure substructure in clusters generally provide estimates of the 
statistical significance of the existence of substructure. The KMM algorithm (Ashman, Bird, & 
Zepf 1994), for example, detects and quantifies the significance that a given cluster is bimodal. 
Mohr, Fabricant, & Geller (1993) introduced the centroid shift to quantify substructure in X-ray 
images of clusters. Although useful for establishing the existence of substructure, these methods 
do not specifically address the relation of the structure to the dynamical state of the cluster. 

We propose a method to quantitatively classify clusters of different morphologies in direct 
relation to the dynamical state of the cluster as indicated by the gravitational potential. We focus 
on structure that is obvious in projection; e.g., the morphological classes described by Jones & 
Forman (1992) for X-ray clusters. That is, the significance of the substructure is a given, what 
the structure implies for the cluster dynamics is the focus of this paper. The method yields 
dimensionless quantities that are especially suited to statistical analysis of a large cluster sample. 
Future papers will apply this technique to a large sample of X-ray clusters (Buote & Tsai 1995b) 
and to clusters generated by N-body / hydrodynamic simulations. In §2. we describe the method. 
We demonstrate the performance of the method on simple models of X-ray clusters in §3. and on 
simulated observations of these models with the ROSAT PSPC in §4.. We apply the method to 
real ROSAT PSPC images of four Abell clusters in §5.. In §6. we discuss the implications of our 
method. In §7. we present our conclusions. 
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2. Method 


The objective of this paper is not to introduce another technique to detect subtle 
manifestations of substructure in galaxy clusters (e.g., Ashman et al. f994; Bird & Beers f993; for 
a comprehensive review see Bird f993), but instead is to propose a simple quantitative scheme to 
categorize clusters (in projection) of different morphological types; e.g., a smooth single-component 
cluster or a widely separated bimodal in the plane of the sky. We focus on subclustering that 
is easily discernible in projection and dynamically important to the whole cluster; i.e. hidden 
substructure along the line of sight is not our concern nor are individual galaxies and small groups 
that do not significantly contribute to the potential energy of the cluster. More precisely, we will 
consider a cluster having two or more components that constitute a sizeable fraction f0%) of 
the total cluster gravitational potential within ~ f Mpc of the cluster center (in projection) to 
possess the type of substructure relevant to this paper; these criteria should include much of the 
substructure suggested by RLT as relevant for cosmological studies. We define a cluster possessing 
this type of substructure to be a multicluster. 

The natural basis for classifying multiclusters is the two-dimensional multipole expansion 
of the projected gravitational potential. Let S(i?,())) be the two-dimensional projection of the 
multicluster mass density, where (i?, <(>) are the conventional polar coordinates. This S(i?,())) 
generates the two-dimensional potential T(i?,())), 

V2T(i?,())) = 47rG'S(i?,())), (f) 


where is the two-dimensional Laplacian and G is the gravitational constant. The standard 
analysis using Green’s functions (e.g., Jackson f975) shows that the potential due to material 
interior to R is, 

- j - 2G ^ {Um cos me/) + bm sin me/)) , (2) 

where the moments and are given by. 


vl/(i?,())) = —2G'aoln 


ajn(R) = / Ti{x') {R')'^ cos, m(j)'(Px', 

Jr’<r 

byn{R) = [ T,[x') (R')"^ sin m(f)'d^x', 

Jr’<r 

and x' = (i?', ()>'). This expansion, similar to its cousin in three-dimensions (e.g., Binney & 
Tremaine f987), has the general properties that (f) the circularly symmetric monopole term (i.e. 
the logarithmic term in eq. 2) is always at least as important as higher-order terms, and (2) the 
dipole (m = f) vanishes when the origin of the coordinate system is set to the center of projected 
mass. T does not represent the total gravitational potential due to S since eq. (2) neglects the 
mass exterior to R. However, as indicated below, we will only be concerned with the gravitational 
effects due to the interior mass. 



A single elliptical cluster only contributes to even terms in the multipole expansion for 'h (eq. 
[2]) when the origin is defined to be the center of projected cluster mass Hence, a significant 
contribution to odd multipole terms unambiguously reflects substructure (asymmetry) in the 
projected cluster, although a multicluster need not have odd multipoles (e.g., a bimodal cluster 
composed of equal-sized subclusters). The even multipoles are also important since multiclusters 
of different morphologies differ in their relative contributions to the even (and odd) multipole 
terms. In keeping with the above definition of a multicluster, we will only consider the first few 
multipole moments {m = 0,1, 2, 3,4) since higher-order terms reflect smaller-scale, dynamically 
less significant structures (see §3.). 

Each term of the multipole expansion (eq. [2]) is a function of position (Rjcj)). Since we want 
to characterize the dynamical state of a large region of a multicluster, a simple procedure is to 
evaluate the multipole moments in a circular aperture; a circular aperture does not introduce 
biases and systematic effects inherent when the aperture shape is modified (through iteration) 
to conform to the shape of an individual multicluster. By computing the moments in a circular 
aperture of radius Rap, the multipole expansion is sensitive to structures having a scale ^ Rap 
and is most sensitive to scales ~ Rap', e.g., if Rap is much greater than any scale associated with 
the cluster then the only significant term in the multipole expansion will be that corresponding 
to m = 0. Since we are uninterested in structure on scales greater than the aperture size, it is 
sensible to neglect the contribution to the potential of mass exterior to Rap as we have done in the 
multipole expansion of eq. (2). 

Because the m > I terms of eq. (2) vanish when integrated over cf), we instead consider 
the magnitude of the terms of each order integrated over cf). Let Tm equal the mth term in the 
multipole expansion of T. Define the quantity 

i,^ap) — X / i,^ap, 0) i,^ap, 0) d()>. (3) 

27r Jo 

Only terms for which m = m' are non-vanishing. Therefore = Pm,m measures the “power” 
within Rap of the terms of order m. Ignoring factors of 2G, these are given by, 

Po = [ao In (Rap)]^ , (4) 


for m = 0 and 

for m > 0. The total power of all multipole moments is simply given by 

OO 

P=Y.Prn. 

m=0 


( 5 ) 

( 6 ) 


^We will, however, make some use of the dipole term (see the following sections). 
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The ratio Pm/P reflects the contribution of the mth multipole moment to the power of the 
total gravitational potential within Ra-p- Since the gravitational potential is directly related to 
the dynamical state of a cluster, the Pm are precisely the measures for substructure we seek. 
Each Pm has units [mass)^ and the scale is set by Rap- Since we are primarily interested in the 
dimensionless ratios of Pm (see §3.) the only dimension that needs to be specified is the scale Rap- 

Because we ignore structures outside of the aperture of radius Rap we naturally may obtain 
different indications for the dynamical state of a cluster depending on the aperture size. For 
example, a widely separated bimodal of equal-sized components may appear essentially relaxed if 
the aperture is small enough so that it encloses only one of the subclusters. If the aperture is large 
enough to include both components, then the Pm will show that the cluster is unrelaxed (actually 
Pm/Po ^ see §§3. and 5.). The ability to quantify the dynamical state of the cluster on varying 
scales is a great asset of the method. If a cluster is virialized then its dynamical state is trivially 
well defined; i.e. the cluster is relaxed on scales larger than the constituent galaxies. However, 
for clusters with significant substructure (e.g., the above bimodal with individually virialized 
subclusters) the dynamical state is a meaningful concept only when referred to a particular scale; 
the relevant scale for the above bimodal is determined by the relative separation of the subclusters. 
The cluster on scales smaller than this, or significantly greater than this, may well be virialized. 

If we know the projected cluster mass density S, then the physical interpretation of the 
Pm is manifest. One may non-parametrically construct an accurate map of S by analyzing the 
weak distortions of background galaxies (e.g., Kaiser & Squires 1993). Unfortunately, since this 
technique is in its infancy and the measurement is difficult (requires sub-arcsecond seeing), such 
maps of the projected cluster mass density exist for only a few clusters (e.g., Fahlman et al 1994; 
Small et al. 1995). In any event, the required source-lens-observer distances (0.15 ^ z ^ 0.6, 
where 2 : is the redshift) fundamentally limits the number of clusters for which weak-lensing maps 
of S may be obtained. 

In order to analyze the structures of a large sample of clusters we must appeal to either 
X-ray images or galaxy positions as a practical substitute for S. X-ray maps of clusters have 
several advantages over galaxy positions. First, the statistical uncertainties associated with galaxy 
positions are intrinsically fixed by the Unite number of cluster galaxies whereas the noise of X-ray 
images is limited only by the sensitivity of the detector and the exposure time of the observation. 
Second, the projection along the line of sight of small groups of galaxies not associated with the 
cluster are less important for X-ray images because the ratio of X-ray luminosity of a cluster to a 
group is much greater than the corresponding ratio of projected galaxy number densities. Third, 
in a similar manner, the projection of small groups within the cluster itself is less significant 
in the X-rays than in the galaxy counts; i.e. the large clumps that are dynamically important 
to the cluster as a whole are proportionally more important in the X-rays. Although X-ray 
images are more suitable for our categorization of multiclusters than galaxy positions, neither 
exactly represents the projected mass density of the multicluster. As a result, we must revise our 
interpretation of the multipole expansion (eq. [2]) when analyzing X-ray images. 
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The X-ray surface brightness, physically differs from S of clusters in several respects. 
Since the X-ray gas emissivity is proportional to p^g^s-, small increases in the gas density enhance 
the X-ray emission correspondingly more than in the underlying mass density. Moreover, the 
relationship of the gas distribution to the underlying mass is unclear because the gas may trace 
the underlying matter, the underlying potential, or neither (e.g., Buote & Tsai f995a). In any 
event, the gas is unlikely to dominate the cluster gravitational potential. 

To achieve successful quantitative classification of multiclusters we require that multiclusters 
of different morphological types as determined from analysis of the projected mass density (S) 
be categorized into corresponding morphological types when using X-ray images. However, the 
actual values of the derived from the X-ray images need not be the same as those derived 
from S. In order to investigate this issue let us consider the projected gas emissivity, as the 
source term of a hypothetical two-dimensional X-ray emissivity potential Consider a single 
elliptical cluster and a widely separated bimodal multicluster with different sized components. 
The qualitative features of the multipole expansion for and T are the same: the odd multipole 
terms will vanish for the single cluster but not for the bimodal multicluster and the relative 
proportions of the will differ for both clusters. Therefore, if the X-ray emission for a given 
multicluster exhibits the same qualitative structure present in the projected mass density S, then 
the multipole analysis of X-ray images should enable quantitative classification of multiclusters 
in the same manner as multipole analysis of S. We explore the feasibility of this scheme in the 
following sections. 


3. Models 


We investigate the capacity of multipole decomposition of X-ray images of galaxy clusters to 
differentiate multiclusters from single-component clusters; by “single-component” we mean one 
cluster component dominates the gravitational potential in projection. To this end we examine 
simple models that capture the features of real X-ray clusters essential for evaluating the utility 
of multipole analysis. Here we shall only consider the fundamental ability of multipole analysis to 
distinguish the different models. In the next section we address the expected real performance 
of the technique by analyzing simulated observations of these models using the instrument 
parameters of the ROSAT PSPC. 

X-ray images of clusters exhibit a variety of morphologies (see Forman & Jones 1990; Jones 
& Forman 1992) ranging from smooth, single-component clusters with regular, nearly elliptical 
isophotes to multi-component clusters possessing several independent emission peaks. We will 
focus our attention on the simplest multicluster, the bimodal. A multicluster with more than two 
components typically has more power in higher order moments and thus is easier to distinguish 
from a single-component cluster; we will, however, investigate real clusters spanning the observed 
range of X-ray morphologies in §5.. Therefore, by examining bimodal multiclusters we will obtain 
a conservative estimate of the viability of multipole analysis. 
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We employ simple models of the aggregate structure of X-ray clusters. For our purposes 
the radial surface brightness profiles of X-ray clusters are sufficiently well parametrized by the 
/3-model (e.g., Cavaliere & Fusco-Femiano f976; Jones & Forman f984; Sarazin f986), 

/ x2]-3d-tl/2 

S^(r) cx f -h [ —] , (7) 

where is the core radius and /3 the slope parameter. From analysis of a sample of bright 
Einstein clusters, Jones & Forman found /3 ~ 0.5 — f and ~ 50 — 750 kpc. Jones & Forman 
noticed that the smooth, single-component clusters with a dominant central galaxy had the 
smallest core radii in their sample, ^ 300. The clusters with large core radii probably possess 
significant core substructure like Abell 2256 (see Jones & Forman f99f; Briel et al. f99f). As 
a result, we assign parameters consistent with the smooth single-component clusters having a 
centrally dominant galaxy to each subclump of a bimodal multicluster. In order to incorporate 
models having a constant ellipticity (e) and orientation, we substitute for r the elliptical radius a, 
where J- jq^, where q is the constant axial ratio. 

We compute on a large sample of cluster models. The single-component clusters are 
represented by single /3-models with e = O.I — 0.6, = fOO — 700 kpc, and /3 = 0.5 — I. 

The ellipticity range spans plausible values for a single-component, ellipsoidal, non-rotating, 
self-gravitating mass (Merritt & Stiavelli 1990; Merritt & Hernquist 1991) appropriate to the case 
where the gas traces the underlying mass. If the X-ray gas is in hydrostatic equilibrium, then 
it will trace the potential and thus be rounder e = 0 ~ 0.3; for a discussion of the evolution of 
X-ray gas shapes see Buote & Tsai (1995a). Although large and e probably characterize a 
multicluster, we include these values for single-component clusters so as to make the distinction 
between multiclusters and single-component clusters more difficult to observe; in this manner 
we may obtain a conservative evaluation of the method. For the bimodals we consider for each 
component = fOO — 300 kpc and for convenience set each /3 to 0.75. The relative separation 
of the components ranges from 250 kpc to 1.5 Mpc and their relative normalization ranges from 
I:I to 1:100. We also allowed the primary component in the bimodal to have e = 0 — 0.4 while 
the secondary component is always kept circular and placed on the major axis of the primary 
component; note all models were normalized to have the same value of uq within an aperture of 
radius 2 Mpc. The range of models we have selected should bracket most bimodals that would be 
considered multiclusters. 

For each model we compute the powers, (eq. [4] and [5]), for m = 0 — 6 within a circular 
aperture of radius Rap = f Mpc; we also examine Rap = 2 Mpc. The projected gravitational 
potential given by equation (2) is defined up to a constant which we specify by choosing units for 
Rap', note the choice of units for Rap is irrelevant for the m > I terms. In order to ensure that the 
Pm have the same values for a cluster independent of its distance we express Rap in units of kpc. 
In Figure I we display contours of four models placed at 2 : = O.IO (see §4.). 

Each choice for the location of the aperture center gives rise to a different multipole expansion 
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(eq. [2]). First, we place the origin of the aperture at the centroid; i.e. where Pi vanishes. By 
doing so, any information possibly found in the first moment is transferred to the higher order 
terms. The quadrupole power, P 2 , in this case, is related to the ellipticity of the cluster. The 
next moment, P 3 , is sensitive to bimodal structure where the components are of unequal size. 

P 4 is similar to P 2 but is more sensitive to smaller-scale structure. Hence P 2 , P 3 , and P 4 yield 
complementary information on the structure of the cluster. 

We investigate a second case where we place the center of the aperture at the peak of the 
surface brightness. To avoid confusion with the previous case, moments generated with the 
origin at the emission peak are denoted Here, the first non-symmetric moment 

only vanishes if the cluster exhibits reflection symmetry about two orthogonal axes centered on 
the origin (such as for a pure elliptical cluster). This moment is then particularly sensitive to 

(pk) 

bimodal multiclusters having nearly equal-sized components. P| ' essentially characterizes a 
circularly-averaged, centroid-shifting power of the cluster within Rap- Although there exists some 
combination of the P^ (centroided) that contains the information given by p[^^\ we employ 
because of its higher sensitivity to nearly equal-sized bimodal multiclusters relative to P 2 , P 3 , and 
P 4 . 

We list in Table f the powers P2, P3, and P4 computed for a selection of the above models 
expressed as a ratio of Pq. These values are indicated by “True” in the table. In Table 2 we list 
the power P^^'^ in terms of Pq^^\ where Pq^^^ is the m = 0 power computed in the circle of 
radius Rap centered at the emission peak. We prefer to consider power ratios, P^/Po, instead of 
the individual P^, because (I) dividing by Pq normalizes to the flux within Rap which enables 
consistent comparison between clusters of different X-ray brightnesses, and (2) the ratio more 
easily distinguishes between an image that is centrally concentrated (i.e. lower Pm/Po) to one 
that is more spread out (i.e. higher Pm/Po); e.g., even if a single and bimodal cluster (of the 
same luminosity) have the same P2 within Rap, the bimodal will necessarily have a smaller Pq. 
Moreover, clusters that are more relaxed should be more dominated by their monopole terms. We 
do not list results for the higher order moments {m > 5) since they are, as expected, particularly 
sensitive to the bimodal multiclusters having relative normalizations greater than 10 : I. These 
models are inconsistent with our definition of a multicluster or the substructure described by RLT 
and are also more susceptible to the hazards of real data (e.g., noise, unresolved sources - see next 
section). 

As expected the odd terms yield the largest differences between single and bimodal clusters 
since they vanish identically for the single-component clusters. The results also demonstrate 
the complementary behavior of j and P3/P0; i.e. the former is larger for equal-sized 
bimodals while the latter vanishes; the former is much less sensitive for unequal-sized bimodals 
while the latter is most sensitive to them. The even ratios, P2/P0 and P4/P0, also perform well at 
distinguishing bimodals from low-ellipticity (e ^ 0.3) single-component clusters. However, there is 
considerable overlap of the flatter single-component clusters (e 0.3) and moderately separated 
bimodals (0.5 Mpc separation of each subclump). For I Mpc separation the degeneracy is lifted 



9 


and the even ratios separate well the single-component and bimodal models. 

From the results listed in Tables 1 and 2, it is clear that the power ratios {m < 4) 
easily discriminate bimodal multiclusters from single-component clusters; here we considered 
multiclusters to have relative separations 500 < r < 1500 kpc and relative normalizations < 10 : 1. 
Moreover, they do not significantly distinguish bimodals with relative separations < 500 kpc 
or relative normalizations > 10 : 1 from single-component clusters implying that these power 
ratios are essentially only sensitive to structure relevant to multiclusters; although we mention 
that P3/P0 and P4/P0 are sensitive to some models with larger relative normalizations when the 
components are sufficiently well separated. 

The ability of the power ratios to distinguish models is dependent on the aperture size. In 
Table 3 we list the power ratios as a function of aperture size (radii 0.5,1.0,1.5, and 2.0 Mpc) for 
two models of very different morphologies: a single-component cluster with e = 0.3, core radius 
300 kpc, and /3 = 0.75; and a bimodal where both components have e = 0, core radius 300 kpc, 
fj = 0.75, are separated by 1 Mpc, and have relative normalization 2 : 1. The power ratios of 
the single-component model (which naturally has vanishing odd P^) are largest for the smallest 
apertures and decrease monotonically with increasing aperture size. This behavior simply reflects 
the decay of higher-order multipole terms with distance just as if the X-rays represented a true 
mass distribution (see §2.); i.e. the qualitative picture of the multipole expansion is preserved 
for our simple models. In the 0.5 Mpc aperture, the bimodal cluster has smaller values of P2/P0 
and P4/P0 than the single-component cluster because the aperture only encloses the dominant 
subcluster; i.e. these power ratios are similar to low-ellipticity single-component models (see 
Tables I and 2). However, the small non-zero values of the odd power ratios, and 

P3/P0, demonstrate that some indication of a large-scale asymmetry is detected within the 0.5 
Mpc aperture. At I Mpc the centroided power ratios peak and then decay for larger apertures 
illustrating the power-ratios’ sensitivity to the scale of the substructure (i.e. I Mpc separation of 
the two components). As expected, the higher-order moments decay most rapidly. In contrast, 
p(p^) Ip(p^) peaks at 1.5 Mpc but otherwise behaves similarly to the centroided moments. 


4. Simulated Observations 


Although we have demonstrated the ability in principle of the power-ratio method to 
quantitatively differentiate single-component clusters from multiclusters, we have not shown if 
these statistics can by usefully constrained by present observations. We now assess the practical 
feasibility of the power-ratio method by simulating real X-ray observations of the cluster models 
with the ROSAT Position Sensitive Proportional Counter (PSPC); for a description of ROSAT 
see Triimper (1983), of the ROSAT X-ray telescope see Aschenbach (1988), and of the PSPC 
see Pfefiermann et al. (1987). The PSPC is more suited to this kind of study than the Einstein 
Imaging Proportional Counter because of its superior resolution (~ 30^' FWHM, on axis) and 
sensitivity. The sensitivity of the PSPC even outweighs the better spatial resolution (~ 4" FWHM) 
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of the ROSAT High Resolution Imager (HRI). Moreover, many more clusters have been observed 
with the PSPC (see § 6 .) than the HRI. The PSPC is also well suited for our study because the 
energy band pass and spectral response implies that the observed depends almost exclusively 
on and is independent of temperature. Hence variations in will not be due to temperature 
fluctuations in the gas. 

Since our primary motive for introducing the power-ratio method is to facilitate consistent 
comparison of the morphologies of statistical samples of clusters (see § 1 . and § 6 .), we simulated 
observations of a well-defined cluster sample. We also chose our observational parameters 
to take advantage of the large number of cluster observations that are becoming available in 
the ROSAT archive. In keeping with the spirit of the power-ratio method (see beginning of 
§2.), we concentrated on bright (for S/N) and nearby {z ^ 0.2; for resolution) clusters. We 
considered clusters having characteristic X-ray luminosity L* obtained from fitting the Schechter 
luminosity function to Abell and AGO clusters in the ROSAT All-Sky Survey (RASS; Ebeling 
1993). Ebeling obtains L* = 1.69 ± 0.37 X 10'^'^ erg s“^ for his sample which translates to a flux 
F* = 1 . 62 :“^ X 10“^^ erg cm“^ s“^, where 2 is the cluster redshift and we have assumed a linear 
distance-redshift relation; this assumption amounts to assuming a static Euclidean universe which 
is suitable for our purposes to explore the effects of angular scale and S/N on the power ratios 
of low-redshift clusters. We placed T* clusters at a series of redshifts (2 = 0.05,0.10, and 0.20) 
appropriate for most of the bright clusters that will be available in the ROSAT archive; a Hubble 
constant Hq = 80 km s“^ Mpc“^ is assumed so these redshifts correspond to 187.5, 375 and 750 
Mpc respectively. To agree with typical observations, we simulated observations having exposure 
times (texp) of 5ks, lOks, and 20ks. 

We constructed a simulated observation for a given model by first choosing the redshift of 
the cluster and then the exposure time. The flux was converted to PSPC counts in the hard 
band (0.4 - 2.4 keV) using the energy-conversion factors (ECFs) given in NRA 9I-OSSA-3, 
Appendix F, ROSAT mission description. For each cluster we assumed a thermal line spectrum 
for temperatures T = 4 — 8 keV and column densities Nh = 10^®'® — 10^^ cm“^ which translates 
to a typical ECF = 0.575 X 10^^ counts cm^ erg“^; the total counts were then {F/.){ECF){texp)- 
The counts within I Mpc of the centroid were normalized to the above total counts. We added 
a uniform background with a count rate of 3 X 10““^ counts s“^ arcmin“^ which is the average 
of the background rates of two PSPC images previously studied (Buote & Canizares 1994,1995). 
The models were then convolved with the off-axis point spread function of the PSPC evaluated at 
I keV (MPE/OGIP Calibration Memo CAL/ROS/93-OI5); note we did not include the support 
structure of the PSPC and thus we confined our study to within the inner 40' diameter ring of 
the PSPC. Moreover, we did not add any exposure variations or vignetting since we assumed the 
observer can adequately correct for these effects. We set the pixel scale to 15" for all the models. 

We included the effects of noise from point sources and Poisson statistics in the simulations. 
First, to each image we added point sources having spatial properties consistent with the PSF of 
the PSPC and numbers consistent with the logW(> S) — log S' distribution given by Hasinger 
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(1991) from analysis of the RASS; see Soltan & Fabricant (1990) and Mohr, Fabricant, & Geller 
(1993) for the inclusion of point sources in simulated Einstein images. The flux of each source 
was determined randomly from the probability distribution dN {> S)/dS and then positioned 
randomly in the field. As before, we converted fluxes to PSPC counts using the ECFs in NRA 
91-OSSA-3. We modeled the point sources with power-law spectra having indices 1.1 — 1.9 and 
Nff = 10^®'® — 10^^ cm“^ which gave ECFs ranging from (0.40 — 0.55) X 10^^ counts cm^ erg“^. 
For each point source we randomly selected an ECF from these values. Any point source that 
would be detected by an observer we excluded from the image, although we did count such a 
source toward the total number of sources dictated by N[> S). We evaluated the significance of 
each source by comparing the source counts to the total counts in a circle of radius 30^' centered 
on the source. We conservatively excluded any source > 5(7 above the total noise. 

For each model corresponding to a particular redshift and exposure time we generated 
1000 simulated observations, and for each simulated observation we computed the power ratios 
P 2 /-P 0 J P 3 /P 0 , and P 4 jPo as described in §3.; here we set Rap = 1 Mpc. Although 
a uniform background will not contribute to for m > 1, it does contribute to Pq. Hence, we 
subtracted the mean background to ensure proper normalization of the power ratios. To evaluate 
-P 2 /-P 0 J P3/P0, and P4/P0 we computed the centroid within Rap and iterated until the centroid 
shifted by less than 10“® pixels (or 25 iterations had been performed). 

We were compelled to adopt a more sophisticated procedure for j because the 
position of the emission peak is much more sensitive to noise then the centroid. For determination 
of the peak position we first convolved the simulated image with a Gaussian filter having a 
corresponding to a physical dimension of 40 kpc. (For example, this corresponds to 3 pixels for 
2 : = 0.05.) Then we took the highest 25% of the pixels within a circle of radius 0.15 Rap centered 
at the peak of the smoothed image and computed the centroid of the highest points. We defined 
this centroid to be the emission peak. Using this position for the peak, we computed pP^'i jpP^'i 
on the original unsmoothed simulated image. 

We computed simulated observations of a subset of models from §3. possessing the essential 
distinguishing characteristics of the single-component and multicluster models. The single¬ 
component models have e = 0, 0.3, and 0.6 while the core radius was set to 300 kpc and /3 to 0.75. 
The multiclusters have for both cluster components e = 0, core radii set to 300 kpc, /3 set to 0.75, 
relative separations 0.5 and 1 Mpc, and relative normalizations 1:1, 2:1, and 5:1. In Figure 2 we 
display contour plots of simulated observations for four of these cluster models. 

The mean value and the 90% confidence limits of the power ratios computed for the 1000 
simulated observations of each model are listed in Tables 1 and 2; for lack of space we only include 
the texp = lOks exposures. The 90% confidence limits are specified so that 10% of the simulated 
observations give power ratios above the upper limit and 10% give power ratios that are below. It 
is clear that for our chosen sample of bright, nearby clusters the power ratios perform nearly as 
well on the simulated observations as on the exact models (§3.) for distinguishing multiclusters 
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from single-component clusters. In fact, the mean value of P 2 /P 0 is generally within ~ 10% of 
the true values and the 90% confidence limits bracketed the true value in all but a few cases. 
p[^^^ /pjf^\ P 3 /P 0 , and P 4 /P 0 behave similarly, although they are usually only within ~ 25% of 
the true values. In addition, for j and P 3 /P 0 the observations give noticeable non-zero 
values (as high as a few times 10“^ for P 3 /P 0 and a few times I0“® for p^P^'> jp^P^'> ) for some 
models where the true value is identically zero. The point sources and Poisson noise do indeed 
contribute noticeable uncertainty, but the uncertainties are generally less than the systematic 
differences in power ratios between the two classes of models. We mention that the contribution 
to noise in the power ratios is roughly equal between the point sources and Poisson noise. We also 
find that the effects of using pixelized images to determine power ratios are very small for the 
chosen redshifts. We do this by essentially determining the difference between a locally smoothed 
version of the image to the observed pixelized image. 

The results for p^^'^ jp^^'^ listed in Table 2 demonstrate that the systematic effects due 
to locating the emission peak are insignificant with respect to noise, pixelization, etc.. This 
systematic effect is due to noise and pixelization of the image only in the sense that the emission 
peak is difficult to determine accurately from a realistic observation. That is, our procedure for 
determining the center by smoothing and then centroiding on the highest surface brightness pixels 
tends to move the center closer to high surface brightness structures which are near the absolute 
peak. For example, if there are two high emission structures which are not extremely widely 
separated, say in a bimodal cluster, then the center will be moved in the direction of the centroid 
of the image. This leads to a smaller amount of power in the given ratio. This effect does not 
influence cases such as the single elliptical cluster where the center is well separated from other 
emission peaks. We emphasize that given realistic observations, the ratio p^^'^ jp^^'^ is a good 
discriminator between single clusters and bimodals, especially those bimodals with nearly equal 
sized components, and brackets the results of the true models within the ~ 90% confidence level. 

We examine whether the sample of multiclusters may be distinguished from the sample of 
single-component clusters by applying a Kolmogorov-Smirnov test to the individual distributions 
of power ratios. That is, we have 1000 simulated observations of 27 single-component cluster 
models and 54 bimodal multicluster models which form two well-defined cluster samples that 
differ only in their morphology, not in their luminosity or distribution or distances; the bimodal 
models with relative normalizations 2:1 were omitted for p^^'^ jp^^'^ leaving 36 total bimodals 
for that case. We obtained, for example, the distribution of P 3 /P 0 values for each sample by 
selecting P 3 /P 0 randomly from one of the 1000 simulations of each model. We list in Figure 3 the 
cumulative distributions of p[^^^/pjf^\ P 2 /P 01 P 3 /P 0 , and P 4 /P 0 obtained in this manner The 
K-S tests yield probabilities of < 0.001% for p[^^^/P q^^\ 3% for P 2 /P 0 , 0.03% for P 3 /P 0 , and 
0.1% for P 4 /P 0 that the singles and bimodals originate from the same population. Thus the K-S 
test for each power ratio convincingly demonstrate that the two data sets could not have come 
from the same parent population. Although the K-S test is useful for hypothesis testing, in §6. we 
discuss the correlations of the power ratios as a means to provide detailed classification of clusters 
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into different morphological types. 


5. Application to Four Real Clusters 


As a final illustration of the performance of the power-ratio technique, we analyzed ROSAT 
PSPC X-ray images of the four Abell clusters A85, A514, A1750, and A2029. A2029 is a classic 
single-component cluster having regular, moderately flattened (e^, ~ 0.15) elliptical isophotes; it is 
nearby (z = 0.0768) and very X-ray bright (~ ^L*; the fluxes of the clusters are listed in Table 5). 
The other three clusters represent the quintessential examples from the qualitative morphological 
classification scheme of Jones & Forman (1992). A85 (z = 0.0556) has a dominant primary 
component with a small secondary and is one of the brightest clusters (~ |T*). A1750 has two 
components of roughly equal size with a respectable X-ray luminosity (~ §T*). Finally, A514 is 
a complex X-ray cluster as classified by Jones & Forman because it has at least three distinct 
emission peaks and highly irregular isophotes; it is at similar redshift as the others (z = 0.0731) 
but is somewhat fainter in X-rays ( ~ ^L*). Therefore, the four clusters span the range of observed 
X-ray morphologies of clusters and hence serve as convenient benchmarks for demonstrating the 
ability of the power-ratios to quantitatively classify real clusters. 

We prepared the images for analysis using the standard IRAF-PROS software. First we 
searched for spikes in the light curves of the images that signal contamination from solar radiation; 
no statistically significant spikes were found for the images. Next we rebinned the images into 15^' 
pixels corresponding to 512 X 512 fields; i.e the same scale used for the simulated observations in 
§4.. Only data from the hard band (0.5 — 2 keV) were used in order to minimize the blurring 
due to the point spread function (PSF) of the PSPC and the contamination from the X-ray 
background. We corrected for exposure variations and telescopic vignetting by dividing the images 
by the exposure maps provided with the observations. To subtract the background we selected an 
area ~ 35' from the field centers apparently uncontaminated by emission from nearby sources or 
the cluster itself. Typically we identified the background region from examination of the radial 
profile centered on the cluster; we designated the background region where the radial profile 
flattened. This procedure undoubtedly suffers to some extent from contamination by emission 
from the cluster, but the errors are insignificant in relation to the total cluster flux which is all 
that we require (see §4.). In Table 4 we list the observational data for the four clusters; in Figure 
4 we show contour plots of the images. 

The final step in the image reduction is to remove embedded sources from the cluster 
continuum. In keeping with the simulated observations of the previous section, we identified 
and removed “obvious” contaminating point sources from visual examination of the images; note 
we arrived at this subjective procedure because automated techniques like the detect package 
in PROS had difficulty identifying sources located in the continuum of the clusters. Although 
some of these sources that we remove are either noise or intrinsic features of the clusters, 
removing them only serves to smooth out the cluster emission thereby decreasing the power from 
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higher-order moments; i.e. these small effects only make a multicluster appear more similar to a 
single-component cluster. Hence, the quantitative significance of a multicluster as separate from 
a single-component cluster is made more robust. The uncertainties and biases regarding source 
removal are investigated in detail in our companion paper (Buote & Tsai 1995b). In Figure 5 we 
show contour plots of the four clusters with the embedded sources removed. 

In Table 5 we list the power ratios computed on these four images for Eap of 0.5, 1, and 1.5 
Mpc, where we have used Hq = 80 km s“^ Mpc“^ as in the previous section. From examination 
of the simulations in §4. we estimate that ~ 50% uncertainties should reflect the 90% confidence 
limits for the clusters. The clusters are clearly differentiated by their power ratios. For example, 
A514 and A2029 have and P3/P0 values that generally differ by over two orders of 

magnitude and the P2/P0 values differ by one order of magnitude, a discrepancy that is highly 
significant considering the 50% uncertainties. In fact, P^^^/Plf^\ P2/P0, and P3/P0 easily 
distinguish A2029 from the each of the other clusters in each aperture. In the 1 Mpc aperture 
P 4 /P 0 for A514 is one hundred times larger than the value for A2029, although the difference in 
P 4 /P 0 is negligible in the other apertures. 

A85 is a prime example of the type of cluster described by RLT as being relevant for 
cosmological tests (§1.). Moreover, next to A2029, it is the cluster with the least obvious structure 
in our sample. In the 0.5 Mpc aperture the P2/P0 and P4/P0 values for A85 are completely 
consistent with those of the smooth A2029. The other two ratios, however, demonstrate substantial 
discrepancy with A2029. All of the power ratios easily distinguish A85 from A2029 in the 1 
Mpc aperture with the ratio P 4 /P 0 being in marginal agreement. Hence the power ratios clearly 
distinguish A85 from single-component clusters like A2029. 

The example of A85 demonstrates that the power ratios are indeed sensitive to the aperture 
size and classify the clusters according to the scale of the substructure they possess as we 
demonstrated for the toy models in §3. and Table 3. A more extreme case is A1750 where the 0.5 
Mpc aperture only encloses one of the subclusters. The power ratios for A1750 on this scale are 
completely consistent with those of A2029 except for P3/P0. It is because a bridge of emission 
connects the two subclusters that P3/P0 is able to measure an asymmetry signaling the presence 
of the subcluster outside the aperture. Unlike A1750, A514 registers very large power ratios in 
all apertures except for P 4 /P 0 which is only larger in the 1 Mpc aperture. The sensitivity to the 
aperture size in this case is a result of the edge of the 1 Mpc aperture falling right on the two 
subclusters to the West. 

The value of j for = 0.5 Mpc for A2029 in Table 5 appears to significantly differ 
from the simulations of single-component clusters in the previous section. For several reasons, we 
do not take this is as an indication of subtle substructure in the central regions of A2029, although 
such structure is in any case not our primary interest. First, A2029 is known to contain a massive 
cooling flow (Sarazin, O’Connell, & McNamara 1992) which implies that the surface brightness 
of the central few hundred kpc does not follow that of the /3 model used in our simulations. 
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Moreover, Sarazin et al. find complicated structure in the X-ray emission that may be due to any 
of a number of possibilities (e.g., lumpy absorption or magnetic fields) not intimately related to 
the total cluster gravitational potential. Since our simulations were restricted to /3 models, they 
may not have adequately accounted for the variety of surface brightness laws to be found in actual 
clusters, especially in the central regions. Hence the range of values quoted in Table 2 for the single 
clusters may not extend to higher values because of the restricted set of models considered and 
thus the presence of subtle substructure should not be concluded. Despite this slight discrepancy 
with the simulations, the relatively small values of j clearly distinguish A2029 from the 
other real clusters. 


6. D iscussion 

The power-ratio method differs from conventional techniques for analyzing substructure 
in galaxy clusters because it is designed to quantitatively label clusters of different aggregate 
morphologies, not simply to quantify the significance of substructure. Moreover, the power-ratio 
method is motivated by cluster dynamics. If the cluster surface mass density is known (e.g., 
from weak-lensing maps), the power-ratio method classifies structure in direct proportion to its 
contribution to the cluster gravitational potential. For X-ray images of clusters this relationship to 
the gravitational potential is not so clear (§2.), but if the X-rays trace the structure in the surface 
density then the dynamical interpretation is qualitatively preserved; this was demonstrated by the 
success of the power ratios at discriminating between the X-ray cluster models and simulations in 
§3. and §4.. 

We envisage the power-ratio method to quantitatively distinguish clusters where the 
substructure is obvious to the eye. For this purpose the significance of the substructure is a 
given, what this structure implies for the aggregate cluster dynamics is our primary concern. 

For example, a definitive measurement of a non-zero value of P 3 unequivocally demonstrates 
that substructure is significant for the cluster, but only in relation to the other can it be 
determined whether the substructure is meaningful to the cluster on a particular scale (Rap)', the 
same argument applies to and therefore to centroid shifts (Mohr et al. 1993). Thus, methods 
that are particularly suited to locate subtle manifestations of substructure and quantifying its 
significance (e.g., the KMM algorithm described by Ashman et al. 1994; also see Bird 1993) are 
not as well suited as the power-ratio method for quantitatively classifying clusters of different 
morphologies. 

To fully realize the capacity of this technique to distinguish clusters by their dynamics, 
the correlations of the power ratios, rather than each ratio alone, should be analyzed; cf. 
the K-S tests in §4.. We illustrate these correlations from analysis of the projections of the 
vector {P2/Pq, P3/Pq, P4/Pq,) onto the two-dimensional coordinate planes for a sample of 
single-component cluster and multicluster models computed in §3.; we exclude p^P^'> jp^P^'> because 
it necessarily has a trivial correlation with the Pm/Po- The sample of cluster models we now 
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consider is slightly more extensive then shown in Tables 1 and 2 to better illustrate the complete 
range of cluster behavior. We consider single cluster models with e = 0.1, 0.2, 0.3, = 100,300 

kpc, and /3 = 0.75; we include only single-component models with e < 0.3 since flatter X-ray 
clusters probably are not relaxed (e.g., Buote & Tsai 1995a). The bimodals have = 300 kpc, 
fj = 0.75, relative separations 0.5,1 Mpc, relative normalizations 1 : 1,2 : 1,5 : 1, and 10 : 1, and 
e = 0, 0.2 for the primary component. For the bimodal models with an elliptical primary (e = 0.2) 
we consider the cases where the secondary is either aligned along the major axis or the minor axis 
of the primary. In Table 7 we list the true power ratios for these models computed in apertures of 
radii Rap = 1,2 Mpc; again we defined Rap in units of kpc. 

We plot the correlations of the power ratios in Figure 6 for Rap = 1 Mpc; since the singles 
(represented by filled ovals) vanish for P 3 /P 0 they have been placed on the P 2 /P 0 or P 4 /P 0 axes 
to show their range in these ratios. In each case the power ratios display positive correlations 
in the sense that clusters with small relative separations and normalizations inhabit the lower 
left of the plots while the top right is populated by the nearly equal-sized, widely separated 
bimodals; i.e. generally the dynamically “mature” clusters are located near the bottom left of 
the plots while “young” clusters populate the upper right. The tightest correlation is in the 
P 2 /P 0 — Pa/Pq plane for the single-component models where essentially P 2 /P 0 cx P^/Pq. The 
bimodal multiclusters, although following the same general trend, have much larger scatter 
especially when the relative normalization >5:1. The correlations clearly show that the 
multiclusters primarily inhabit a localized region of the three-dimensional space of power ratios; 
this space is typically P 2 /P 0 = 10“® — 10“'^, P 3 /P 0 = f0“® — I0“®, and P^jP^ = I0“® — I0“®. 
Clusters with smaller relative separations 0.5 Mpc) and larger relative normalizations 10 : I) 
generally lie outside this volume in the direction of small values of the power ratios. 

In Figure 6 we draw a dashed box to represent the region inhabited by single-component 
clusters considering the effects of noise, pixelization etc. from simulated observations of these 
clusters with the PSPC (see §4.). Although some of the multiclusters lie in this region in one 
of the planes, they are usually removed from correlation with the remaining power ratios. For 
example, the bimodal with relative separation I Mpc, relative normalization 1:1, and primary 
e = 0.2 with secondary on the major axis is the model with the lowest P 3 /P 0 and P 4 /P 0 shown 
in Figure 6 and within the error box for single-component clusters in the P 2 /P 0 — Pa/Pq plane. 
However, this model lies outside the error boxes in the other planes and is thus distinguished 
as a multicluster. In general, where the clusters lie in the volume classifies them by the type of 
structure they possess; i.e. clusters falling in the multicluster region possess substructure that is 
dynamically relevant to the aggregate cluster dynamics on a scale Rap (see §2.). 

The correlations of the power ratios for Rap = 2 Mpc are displayed in Figure 7. Although the 
power ratios still exhibit positive correlations (with values about 10 times less than before), the 
distinction between the single-component models and the bimodals is not as pronounced as for 
the Rap = f Mpc case - especially upon considering the effects of real observational uncertainty 
(i.e. as above the dotted lines show the allowed region for single-component clusters considering 
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simulated observations in §4.). This is again a manifestation of the multipole description at work. 
That is, since the bimodal models we have constructed have separations either 0.5 or f Mpc, the 
higher-order multipoles are largest on those scales. For aperture sizes larger than these scales the 
monopole term quickly dominates as the higher-order moments rapidly decay. It is thus important 
to examine different aperture sizes to determine on which scale substructure is particularly 
important for real clusters. (We mention that the bimodal model with e = 0.2, Eg = 0.5 Mpc, and 
EEL = 2:1 and secondary aligned along the major axis of the primary has a P3/P0 value larger 
at Eap = 2 Mpc than Eap = f Mpc - although both values are very small. This appears to be an 
interesting case where the ellipticity of the primary acts to reduce the value of P3/P0 which does 
not happen for the models where the secondary is aligned along the minor axis of the primary. 

In a similar case, the bimodal model with e = 0.2, Eg = I.O Mpc, EEL = 10 : I and secondary 
aligned along the minor axis of the primary has P2/P0 value larger at Eap = 2 Mpc than Eap = 1 
Mpc. More interesting, though, is that this model has a smaller value of P2/P0 than P4/P0 for 
Eap = 1 Mpc. We find that these interesting cases are not well-represented by real clusters [see 
Buote & Tsai 1995b] which may simply be the result of the cluster initially collapsing along its 
shortest axis - e.g., Lin, Mestel, & Shu 1965.) 

A potential application of the power-ratio method is for determining quantitatively whether 
a quasi-hydrostatic equilibrium description for the X-ray emitting gas of an individual cluster is 
justified for the purpose of constraining its intrinsic shape and its total mass distribution. Buote & 
Tsai (1995a) tested the viability of X-ray analysis for constraining the intrinsic shapes of clusters 
of galaxies using the simulation of Katz & White (1993). They concluded that at low redshifts 
{z ^ 0.25) the X-ray method accurately measured the true ellipticity of the three-dimensional 
cluster dark matter up to projection effects. At higher redshifts {z ^ 0.25), however, the X-ray 
method yielded unreliable results since the gas does not trace the cluster gravitational potential. 
Buote & Tsai proffer some necessary conditions for the reliability of X-ray methods: (I) that 
there is no obvious substructure on the same scale used to compute the aggregate shape and (2) 
the isophotes are regularly shaped and not too elongated (e^, ^ 0.3). The power-ratio method is 
particularly suited to quantify these necessary conditions. Results on our study will be presented 
elsewhere, but see Buote & Tsai (1995b) for an outline of a prescription for this program. 

The power-ratio method is ideally suited to constrain LI via the Morphology - Cosmology 
connection (see §1.; RLT; Evrard et al. 1993). The method provides a simple, consistent 
comparison of the structure of clusters since the power ratios are (I) computed in a well-defined 
aperture, (2) normalized to the flux within that aperture, and (3) do not require any fitting. 
Because it is particularly sensitive to structure relevant to the dynamical state of the cluster, 
the power-ratio method specifically quantifies the type of substructure described by RLT as 
relevant for cosmology. The ratio P3/P0 (and P4/P0), being sensitive to unequal-sized bimodal 
multiclusters, is most relevant to the structure envisioned by RLT, while j and P2/P0 
are more sensitive to roughly equal-sized subclumps. 

A large number of clusters similar to those of our sample in §4. will be available in the 
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ROSAT archive. From examination of the ROSAT master log of pointed observations (in the 
HEASARC-Legacy database) for Abell clusters having (1) measured flux 10“^^ erg cm“^ s“^ 
as published by Ebeling (1993), (2) exposure times > 5ks, and (3) 2 : ^ 0.2 we find ~ 50 eligible 
clusters. Higher redshift clusters will be available for analysis with AXAF because of its superior 
resolution. As a result, at least 124 Abell clusters from Ebeling (1993) having flux > 10“^^ erg 
cm“^ s“^ will in principle be eligible for analysis. Thus application of the power-ratio method 
to these samples should enable a thorough statistical investigation of the viability of using the 
observed structure of clusters to place interesting constraints on 12. In Buote & Tsai (1995b) we 
apply the power-ratio method to 55 ROSAT PSPC clusters to furnish a catalog of power ratios 
suitable for statistical analysis. 

An effort to analyze the Morphology - Cosmology connection was undertaken by Evrard et 
al. (1993). These authors proposed to quantify the structure of Einstein clusters using a mean 
centroid shift (Mohr et al. 1993), a mean axial ratio, and a mean slope of the surface brightness 

(vk) 

of the entire X-ray image; the mean centroid shift and mean axial ratio are related to our a) , 

(vk) 

5) ' and 02 , 62 - By computing the power ratios in apertures defined by the cluster distances, we 
consistently sample the same intrinsic scales of clusters. In addition, by using a series of aperture 
sizes we also obtain information regarding the scale of the substructure in the cluster sample. 
Evrard et ah, in contrast, sample different cluster scales for each cluster because they compute 
mean quantities for the entire X-ray images; the size of a cluster X-ray image is dependent on the 
flux, intrinsic size, and distance of the cluster. Evrard et al. also do not employ a third (or forth) 
moment which is in fact more sensitive to the unequal-sized bimodal multiclusters envisioned by 
RLT than the lower-order moments. 

The manner in which Evrard et al. compute the centroid-shift and axial ratio (as explained in 
Mohr et al. 1993) also makes a direct comparison to intrinsic properties of the cluster uncertain. 
Specifically, the surface brightness in a circular annulus of a given width is first Fourier expanded 
to low order given a trial center for the annulus. The expansion is then fitted to the image taking 
the Fourier coefficients as parameters of the fit. The location of the center of the image is then 
iterated so as to minimize the coefficient of the m = 1 term (Ci). This center is then used to 
compute the centroid shift and the axial ratio. This procedure gives correct values for these latter 
quantities only if the cluster being considered is very nearly elliptical and higher order terms in the 
Fourier expansion of the surface brightness are small. This is because the values of the coefficients 
determined by fitting a highly truncated version of the expansion to the image are not necessarily 
the true values of the Fourier coefficients when higher order terms are important. Given a center, 
the true Fourier coefficients are given by moments of the surface brightness distribution; fitted 
values will depend on the highest order considered and will not have their usual meanings. 

Mohr et al. (1993) only considered clusters for which the condition that higher order terms in 
the expansion of the surface brightness be small is satisfied. Their values of the centroid shift and 
the axial ratio are probably accurate. However, consideration of substructure in these clusters are 
of limited use in cosmological considerations since the substructure required by RLT will indeed 
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give rise to significant higher order terms, as seen in §3. 


7. Conclusions 


We have described a technique to quantitatively classify clusters of galaxies according to their 
projected morphology. In particular, we addressed structure that is easily discernible in projection 
and dynamically important to the whole cluster. A cluster possessing substructure of this type 
(§2.) we defined to be a multicluster. We specifically designed our method to quantitatively 
distinguish multiclusters from single-component clusters in projection. The method is derived 
from the two-dimensional multipole expansion of the projected cluster gravitational potential; 
i.e. the square of each multipole term averaged over a circular aperture is called the power, P^, 
and when divided by another (particularly Pq) we call it a power ratio. For the case where 
the surface mass density of the cluster is known, e.g., from analyzing the weak distortions of 
background galaxies (e.g., Kaiser & Squires 1993), the power-ratio method classifies structure in 
relation to its contribution to the cluster gravitational potential. 

For X-ray images of clusters this relationship to the gravitational potential is not so 
transparent (§2.), but if the X-rays approximately trace the structure in the surface mass density 
then qualitatively the dynamical interpretation is preserved. We demonstrated this assertion by 
analyzing the performance of the power-ratio method applied to simple models capturing the 
essential features of real X-ray clusters. In particular we focused on models of single-component 
clusters and bimodal multiclusters. By construction the structure of these models of X-ray 
clusters reflected structure in the projected mass; clumps in the X-rays corresponded to the same 
clumps in the mass, although not in exactly the same proportions. We determined that the ratios 
P 2 IPq, Pq/Pq, and P 4 /P 0 performed best for distinguishing between single-component 
clusters and multiclusters; the powers were computed in a circle located at the cluster centroid 
except for j which was centered at the emission peak. The ability of the power ratios 
to differentiate clusters was optimized when the aperture size was of order the separation of the 
clumps of the bimodals. 

We simulated observations of these models using the instrument parameters of the ROSAT 
PSPC for a sample of nearby [z < 0.2) and bright (flux = F* = 1 . 62 :“^ X 10“^® erg cm“^ s“^) 
clusters. The effects of point sources. X-ray background, Poisson noise, and realistic exposure 
times were also incorporated into the simulations. The power ratios perform nearly as well 
on the simulated observations as on the exact models (§3.) for distinguishing multiclusters 
from single-component clusters. The point sources and Poisson noise do contribute noticeable 
uncertainty, but generally less than the systematic differences in power ratios between the two 
classes of models; the contribution to noise in the power ratios is roughly equal between the point 
sources and Poisson noise. Applying a Kolmogorov-Smirnov test to the power-ratios obtained 
from the simulations clearly demonstrates the sample of multiclusters may be distinguished from 
the sample of single-component clusters. 
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We applied the power-ratio method to ROSAT PSPC images of A85, A514, A1750, and 
A2029. These clusters each have very different X-ray morphologies with A2029 being a smooth, 
single-component cluster; A85 being a dominant smooth component with a small secondary; 
A1750 being two components of nearly equal size; and A5f4 being the quintessential complex 
cluster in the qualitative classification scheme of Forman & Jones (f990; Jones & Forman f992); 
A85 and Af750 are also listed by Forman & Jones as the definitive members of their own classes. 
We find that the power ratios easily differentiate the clusters, especially when the aperture size 
is f Mpc {Hq = 80 km s“^ Mpc“^). In a companion paper (Buote & Tsai 1995b) we apply the 
power ratios to a large sample of clusters observed with the PSPC. 

We have discussed the suitability of the power-ratio method to constrain 17 via the Morphology 
- Cosmology connection (see §1.; RLT; Evrard et al. 1993). The method provides a simple, 
consistent comparison of the structure of clusters since the power ratios are (I) computed in a 
well-defined aperture, (2) normalized to the flux within that aperture, and (3) do not require 
any fitting. Moreover, the power-ratio method is specifically sensitive to the type of substructure 
described by RLT as relevant for cosmology. We also discussed the ability of the power ratio 
method to assess the viability of a particular cluster being described by hydrostatic equilibrium 
for the purposes of X-ray analysis of its intrinsic shape and of its total mass distribution. 
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Table 1. Power Ratios 


Model True z = 0.05 


z = 0.10 


z = 0.20 


Ti/A (lO-O 


Single (e): 


0 

0.000 

0.160 

0.068 - 0.182 

0.482 

0.181 - 0.548 

1.61 

0.609 - 1.86 

0.3 

39.6 

41.9 

38.4 - 45.4 

42.5 

35.6 - 49.4 

43.2 

29.9 - 56.8 

0.6 

162 

167 

160 - 175 

169 

155 - 183 

168 

141 - 196 

Bimodal { Rs , 

REL ): 







(0.5,1:1) 

59.9 

62.7 

58.0 - 67.2 

63.0 

53.2 - 72.8 

63.7 

45.0 - 84.0 

(0.5,2: 1) 

46.3 

48.0 

43.7 - 52.1 

49.1 

40.5 - 57.7 

48.7 

32.7 - 66.1 

(0.5,5 : 1) 

16.8 

18.3 

15.7 - 21.0 

18.8 

13.7- 24.0 

19.2 

9.32 - 30.2 

(1.0,1:1) 

882 

865 

845 - 887 

876 

834 - 918 

860 

707 - 945 

(1.0,2: 1) 

550 

525 

507 - 544 

554 

516 - 596 

534 

461 - 614 

(1.0,5: 1) 

85.9 

82.2 

75.1 - 89.1 

90.9 

77.1 - 107 

90.5 

58.9 - 122 


O 

1—1 

-0 







Single (e): 








0 

0.000 

0.030 

0.003 - 0.067 

0.116 

0.012 - 0.275 

0.456 

0.049 - 1.06 

0.3 

0.000 

0.025 

0.003 - 0.058 

0.103 

0.009 - 0.244 

0.413 

0.039 - 0.945 

0.6 

0.000 

0.021 

0.003 - 0.048 

0.081 

0.009 - 0.180 

0.337 

0.031 - 0.767 

Bimodal [ Rg , REL ): 







(0.5,1:1) 

0.000 

0.032 

0.004 - 0.075 

0.133 

0.013 - 0.320 

0.501 

0.054 - 1.14 

(0.5,2: 1) 

0.471 

0.505 

0.291 - 0.727 

0.626 

0.179 - 1.16 

0.985 

0.134 - 2.17 

(0.5,5 : 1) 

0.664 

0.693 

0.441 - 0.960 

0.805 

0.299 - 1.38 

1.18 

0.196 - 2.47 

(1.0,1:1) 

0.000 

0.061 

0.007 - 0.146 

0.275 

0.026 - 0.657 

1.045 

0.112 - 2.44 

(1.0,2: 1) 

29.6 

28.4 

26.2 - 30.6 

28.7 

24.2 - 33.2 

28.6 

20.4 - 37.3 

(1.0,5: 1) 

14.1 

12.6 

11.4 - 14.0 

14.2 

11.4- 17.1 

14.3 

8.54 - 20.2 








Table 1—Continued 


Model True z = 0.05 z = 0.10 z = 0.20 


Pi/Po (lO-O 


Single (e): 


0 


0.000 

0.014 

0.001 

- 0.033 

0.054 

0.006 

- 0.121 

0.209 

0.022 

- 0.474 

0.3 


0.099 

0.135 

0.068 

- 0.205 

0.170 

0.042 

- 0.324 

0.320 

0.046 

- 0.703 

0.6 


2.18 

2.24 

1.97 - 

2.51 

2.34 

1.76 - 

2.95 

2.53 

1.43 - 

3.76 

Bimoda 

{Rs 

,REL): 










(0.5,1 

1) 

0.036 

0.068 

0.021 

- 0.120 

0.107 

0.014 

- 0.233 

0.303 

0.037 

- 0.700 

(0.5,2 

1) 

0.051 

0.080 

0.029 

- 0.139 

0.125 

0.015 

- 0.271 

0.279 

0.031 

- 0.628 

(0.5,5 

1) 

0.056 

0.087 

0.031 

- 0.150 

0.131 

0.020 

- 0.281 

0.273 

0.032 

- 0.610 

(1.0,1 

1) 

7.90 

7.75 

6.78 - 

8.75 

8.12 

6.57- 

9.67 

8.27 

5.37 - 

11.3 

(1.0,2 

1) 

10.1 

9.65 

8.84 - 

10.4 

10.1 

8.36 - 

11.7 

9.94 

6.67 - 

13.3 

(1.0,5 

1) 

3.98 

3.68 

3.23 - 

4.14 

4.15 

3.19 - 

5.16 

4.37 

2.23 - 

6.63 


Note. — The power ratios are computed in a 1 Mpc circular aperture about the centroid. “True” 
corresponds to the intrinsic power ratio from §3. and the values for the different redshifts are the mean and 
90% conhdence limits for the 1000 simulated observations of the models (see §4.) having exposure time lOks. 
The core radii and [3 parameters of the models are hxed as described in §4.. The single-component models 
only differ in their ellipticity e. The bimodal multiclusters models are listed for different values of their 
relative separation, Rg (Mpc), and their relative normalization, REL. 








Table 2. Power Ratios 


Model True z = 0.05 z = 0.10 z = 0.20 


pm/pipk) 


Single (e): 


0 

0.000 

2.36 

0.078 - 3.31 

5.21 

0.370 - 8.43 

14.8 

1.64 - 29.6 

0.3 

0.000 

0.798 

0.065 - 1.49 

2.87 

0.239 - 6.68 

11.6 

1.10 - 25.9 

0.6 

0.000 

0.588 

0.051 - 1.27 

2.47 

0.207 - 3.77 

8.17 

0.758 - 19.3 

Bimodal {Rg, 

REL): 







(0.5,1:1) 

4437 

4378 

3762 - 4851 

3705 

2887 - 4706 

2896 

1336 - 4192 

(0.5,5 : 1) 

472 

444 

294 - 621 

350 

214 - 614 

254 

39.9 - 666 

(1.0,1:1) 

5148 

7092 

6825 - 7318 

5991 

5574 - 6408 

6149 

5289 - 6955 

(1.0,5: 1) 

372 

399 

294 - 512 

322 

202 - 509 

257 

69.4 - 575 


Note. — The power ratio is computed assuming an aperture of radius 1 Mpc centered on the emission 
peak. Quantities are listed as in Table 1 
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Table 3. Power Ratios vs. Aperture Size 


^ap 

(Mpc) 

pipi^) j p^P^) 

P 2 IP 0 

P 3 /P 0 

P 4 IP 0 


Single: e = 0.3: 

0.5 

0. 

93.6 

0. 

0.235 

1.0 

0. 

39.6 

0. 

0.099 

1.5 

0. 

17.8 

0. 

0.038 

2.0 

0. 

9.26 

0. 

0.017 

Bimodal [Rg = 1.0, REL = 2:1): 

0.5 

16.5 

0.822 

0.062 

0.007 

1.0 

1824. 

550. 

29.6 

10.1 

1.5 

3288. 

166. 

3.38 

0.720 

2.0 

2078. 

55. 

0.643 

0.083 


Note. — The power ratios are expressed in units of 10 ^ 


Table 4: Observational Parameters 


Cluster 

z 

Exposure (ks) 

0.1-2.4 keV Flux 
(10“^^ erg cm“^ s“^) 

Background 

(10“'^ cts s“^ arcmin“^) 

A85 

0.0556 

10.240 

80.61 

3.22 

A514 

0.0731 

18.111 

5.00 

2.24 

A1750 

0.0855 

13.148 

14.62 

2.79 

A2029 

0.0768 

12.550 

66.67 

5.50 


Note. — Only the 0.5 and 1 Mpc values are listed for A85 because 1.5 Mpc lies outside the central ring of the 
PSPC. The fluxes are from Ebeling (1993) except A514 which we computed in this paper (see §5.). The background 
rate is computed in regions ~ 30 — 40* from the field centers. 
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Table 5: Power Ratios of Abell Clusters 


Cluster 

pVkjppk 

0.5 1.0 1.5 

0.5 

T 2 /T 0 

1.0 

1.5 

0.5 

T 3 /T 0 

1.0 

1.5 

0.5 

T 4 /T 0 

1.0 

1.5 

A85 

311 

333 


15.4 

13.8 


1.10 

0.811 


0.032 

0.166 


A514 

9369 

4833 

2557 

273 

300 

179 

22.7 

5.83 

7.40 

0.679 

13.6 

0.155 

A1750 

7 

3670 

3996 

8.9 

818 

311 

8.60 

6.08 

0.679 

0.058 

12.6 

7.33 

A2029 

37 

9 

5 

14.0 

1.7 

2.0 

0.031 

0.004 

0.020 

0.050 

0.073 

0.060 


Note. 


Power ratios in units of 10 ^ 


for real PSPC images of Abell clusters computed for aperture radii 0.5, 


1 , 


and 1.5 Mpc. 







Table 6. 


Correlations of the Power Ratios 



1 p(p>^) 

P2/P0 

Ps/Po 

P4/P0 

Models 

f.O 2.0 1 

1.0 2.0 1 

1.0 2.0 1 

1.0 2.0 


Single (e,ax): 

(0.1,0.1) 

0. 

0. 

0.527 

0.081 

0. 

0. 

7.8e-5 

9e-6 

(0.1, 0.3) 

0. 

0. 

4.15 

1.05 

0. 

0. 

9.6e-4 

1.9e-4 

(0.2, 0.1) 

0. 

0. 

2.07 

0.314 

0. 

0. 

1.3e-3 

1.5e-4 

(0.2, 0.3) 

0. 

0. 

17.1 

4.16 

0. 

0. 

0.017 

3.2e-3 

(0.3, 0.1) 

0. 

0. 

4.55 

0.680 

0. 

0. 

7.1e-3 

8.0e-4 

(0.3, 0.3) 

0. 

0. 

39.6 

9.26 

0. 

0. 

0.099 

0.017 

Bimodal [e, Rs, REL): 

(0.0, 0.5,1 : 1) 

4437 

1223 

59.9 

4.44 

0. 

0. 

0.036 

2.3e-4 

(0.0, 0.5, 2 : 1) 

1929 

543 

46.3 

3.50 

0.471 

0.010 

0.051 

3.3e-4 

(0.0, 0.5, 5 : 1) 

472 

135 

16.8 

1.36 

0.664 

0.016 

0.056 

3.9e-4 

(0.0, 0.5,10 : 1) 

139 

40.3 

5.68 

0.480 

0.330 

8.3e-3 

0.030 

2.3e-4 

(0.0, 1.0, 1 : 1) 

5148 

4716 

882 

70.6 

0. 

0. 

7.90 

0.059 

(0.0, 1.0, 2 : 1) 

1824 

2078 

550 

55.3 

29.6 

0.643 

10.1 

0.089 

(0.0, 1.0, 5 : 1) 

372 

515 

85.9 

21.1 

14.1 

0.968 

3.98 

0.096 

(0.0, 1.0, 10 : 1) 

101 

153 

17.5 

7.31 

3.47 

0.502 

0.993 

0.055 

Major Axis: 

(0.2, 0.5, 1 : 1) 

4367 

1217 

95.2 

9.74 

0.480 

0.041 

0.203 

5.7e-3 

(0.2, 0.5, 2 : 1) 

1889 

539 

90.6 

10.4 

3.7e-3 

6.4e-3 

0.184 

5.6e-3 

(0.2, 0.5, 5 : 1) 

460 

134 

56.7 

8.19 

0.179 

1.5e-4 

0.144 

5.3e-3 

(0.2, 0.5,10 : 1) 

135 

39.9 

37.6 

6.48 

0.117 

5.6e-4 

0.092 

4.7e-3 

(0.2,1.0,1 : 1) 

5042 

4690 

993 

88.2 

1.15 

0.152 

12.5 

0.155 

(0.2,1.0, 2 : 1) 

1779 

2062 

677 

76.9 

19.1 

0.202 

13.0 

0.159 

(0.2, 1.0,5 : 1) 

362 

510 

163 

39.3 

10.3 

0.577 

4.55 

0.144 

(0.2, 1.0, 10 : 1) 

98.6 

151 

64.6 

20.7 

2.57 

0.330 

1.25 

0.086 

Minor Axis: 

(0.2, 0.5, 1 : 1) 

4367 

1217 

36.4 

1.25 

0.629 

0.046 

1.9e-3 

7.6e-5 

(0.2, 0.5, 2 : 1) 

1889 

539 

17.4 

0.276 

1.97 

0.086 

0.039 

1.3e-3 

(0.2, 0.5, 5 : 1) 

460 

134 

0.401 

0.285 

1.53 

0.059 

0.099 

3.7e-3 

(0.2, 0.5,10 : 1) 

135 

39.9 

2.01 

1.35 

0.660 

0.026 

0.080 

4.2e-3 

(0.2,1.0,1 : 1) 

5042 

4690 

850 

56.6 

3.19 

0.186 

5.49 

0.020 

(0.2,1.0, 2 : 1) 

1779 

2062 

423 

37.7 

47.1 

1.43 

9.12 

0.062 

(0.2,1.0,5 : 1) 

362 

510 

28.6 

8.35 

17.6 

1.49 

4.09 

0.112 

(0.2,1.0,10 : 1) 

98.6 

151 

0.026 

0.708 

4.29 

0.716 

1.19 

0.078 


Note. — The power ratios (in units of 10~^) computed in apertures of radii f Mpc and 2 Mpc; units of and 
Rs are also Mpc and e in the bimodal models refers to the ellipticity of the primary component. Bimodal models 
where the secondary component lies along the major axis of the primary are listed under the heading “Major Axis”; 
similarly the models where the secondary is aligned along the minor axis are listed under “Minor Axis”. 
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Fig. 1.— 

Contour plots of four simple models for X-ray clusters placed at 2 : = O.fO: (a) is a 
single-component model with core radius 300 kpc and e = 0.30; (b) - (d) are bimodals each 
separated by f Mpc and have core radii 300 kpc but the clumps are in proportion f : f for (b), 

2 : f for (c), and 5 : f for (d). The contours are separated by factors of two in surface brightness 
(arbitrary units). The units are in f5" pixels which translates to 36.7 pixels/Mpc and 2.7 
Mpc/side for Hq = 80 km/s/Mpc. 


Fig. 2.— 

Contour plots of fOks simulated observations of the cluster models in Figure f as described in §4.. 
The images have been smoothed with the PSPC PSF for viewing purposes only. 


Fig. 3.— 

Cumulative distributions of P 2 /P 01 P 3 /P 0 , and P 4 /P 0 for the simulated cluster 

samples of single-component clusters (solid) and bimodal clusters (dotted). 


Fig. 4.— 

Contour plots of the PSPC images of Abell clusters A85, A5f4, Af750, and A2029 corrected for 
exposure, vignetting, and background; the angular sizes of the fields are the same as Figure f. 
The contours are separated by factors of 2 in intensity and the images have been smoothed with 
the PSPC PSF for viewing purposes only. 


Fig. 5.— 

Contour plots of the PSPC images of Abell clusters A85, A5f4, Af750, and A2029 prepared as in 
Figure 4 but with sources removed as described in §5.. The images have been smoothed with the 
PSPC PSF for viewing purposes only. 


Fig. 6.— 

Correlations of the power ratios for a sample of models of single-component clusters (filled ovals) 
and bimodal multiclusters (crosses); see Table 3 for a description of the sample. The dashed lines 
represent the region where single-component clusters are allowed due to observational uncertainty 
as determined by the simulated PSPC observations in §4.. 


Fig. 7.— 

Same as Figure 6 for the power ratios computed in the 2 Mpc aperture. 



